Non-Gaussian Halo Bias Beyond the Squeezed Limit 
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Primordial non-Gaussianity, in particular the coupling of modes with widely different wavelengths, 
can have a strong impact on the large-scale clustering of tracers through a scale-dependent bias 
with respect to matter. We demonstrate that the standard derivation of this non-Gaussian scale- 
dependent bias is in general valid only in the extreme squeezed limit of the primordial bispectrum, 
i.e. for clustering over very large scales. We further show how the treatment can be generalized 
to describe the scale-dependent bias on smaller scales, without making any assumptions on the 
nature of tracers apart from a dependence on the small-scale fluctuations within a finite region. If 
the leading scale-dependent bias Ab oc k", then the first subleading term will scale as fc"+^. This 
correction typically becomes relevant as one considers clustering over scales k > 10^^ hMpc~^ . 
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I. INTRODUCTION 

Primordial non-Gaussianity is one of the most promis- 
ing probes of the physics and nature of inflation in the 
early Universe [1]. Currently the best constraints on non- 
Gaussianity come from the cosmic microwave background 
as observed by the Planck satellite [2]. However, it has 
become clear recently that observations of the clustering 
of large-scale structure (LSS) tracers will offer competi- 
tive constraining power on non-Gaussianity [3~9]. 

The key ingredient in describing the impact of non- 
Gaussianity on LSS statistics is the description of the bi- 
asing of tracers, that is, the relation to the matter density 
field. Dalai et al. [10] showed that in the presence of non- 
Gaussianity of the local type, where the non-Gaussian 
potential perturbation (f> is given in terms of a Gaussian 
field (h via 



(x)-(^(x) + /nl <^'(x)-(0') 



(1) 



where /nl is a dimensionless parameter, leads to a 
strongly scale-dependent bias which increases towards 
large scales as k^^ in Fourier space. That is, on large 
scales, tracers follow the potential, rather than the mat- 
ter as for Gaussian initial conditions. This effect can be 
understood as follows. Given a potential described by 
Eq. (1), one can easily show [10, 11] that the power spec- 
trum of small-scale matter fluctuations in a patch of size 
Rl around x is rescaled as 



P(fc,) = [l + 4/NL0L(x)]P(fc,)^ 



(2) 



where <J)l{'^) is the potential averaged over the patch, 
and P{ka) is the matter power spectrum derived from the 
Gaussian field (j). Tracers in this patch then effectively 



•Einstein fellow 



form in a Universe with higher primordial power spec- 
trum amplitude, which will clearly change their abun- 
dance relative to other patches resulting in a modulation 
of the tracer density by long wavelength potential per- 
turbations. 

The rescaling Eq. (2) is only valid in the large-scale 
limit Rl — > oo, that is when gradients of can be ne- 
glected. As we will see, this is formally equivalent to 
only considering the leading contribution to the bispec- 
trum Brj, in the squeezed limit, lim/j_5.o B^{k, kg, [kg +k|). 
This is clearly not a good assumption for most current 
large-scale structure surveys, which probe Fourier modes 
k > 10^^ /iMpc^^. Several treatments have gone beyond 
the approximation Eq. (2) [12-14], however they all as- 
sumed specific models for the tracers (thresholding or 
excursion set). It is likely that more refined models for 
tracers will be necessary in order to adequately describe 
the large samples of galaxies, quasars, and other tracers 
delivered by ongoing and upcoming surveys. The goal of 
this paper is to go beyond the limit described by Eq. (2) 
while keeping the treatment fully independent of detailed 
assumptions about the tracers. 

Our treatment will be based on the approach devel- 
oped in Schmidt et al. [15] (see also [16, 17]). The under- 
lying idea is that, when coarse-grained on a sufficiently 
large scale Rl , the abundance of tracers only depends lo- 
cally on various coarse-grained properties of the density 
field. The properties considered in [15] were the coarse- 
grained matter density p^ =p{1 + Sl)j the curvature of 
the matter density V^pi, and the amplitude of small- 
scale fluctuations j/* . The assumption of locality is valid 
as long as Rj^ is much larger than the scale of non-locality 
of the tracer considered. By defining renormalized bias 
parameters, it is then possible to absorb the dependence 
on the arbitrary scale Rl in the expression for tracer cor- 
relations, which then becomes manifestly convergent in 
terms of large-scale matter correlators. The renormalized 
bias parameters are then given in terms of derivatives of 
the mean abundance of tracers with respect to the prop- 
erties of the background Universe and initial conditions. 



Let us briefly recap the renormalization procedure pre- 
sented in [15] for tracer clustering in the presence of 
primordial non-Gaussianity. In this paper we will as- 
sume that the leading non-Gaussian contribution is given 
by the bispectrum; the generalization to higher A'^-point 
functions is straightforward. If we assume that, on large 
scales, tracer correlations are completely described by the 
dependence of the tracer density on the coarse-grained 
fractional matter density perturbation Sl, 



-X^)= j d?yWL(\^-y\)5{y), 



(3) 



where Wl is an arbitrary spherically symmetric filter 
function on the scale Rl , then the tracer correlation func- 
tion ^Ti is to leading order given by 



Ur)^hlUr) + bib2{5L{l)5l{2)) 



(4) 



where S,L{f) ~ ('5L(l)(5i(2)), and h^ are the renormal- 
ized "PBS" bias parameters defined as derivatives of the 
mean tracer abundance with respect to a change in the 
background density p of the Universe. 

The leading non-Gaussian modification is the sec- 
ond term oc &i62- In terms of the matter bispectrum 
Bm{ki, ^2,^3), this three-point correlator is given by 



{5l{1)61{2)) = 



(Pk 



Wdk) 



(27r)3 
X WL{ki)WL{\k 



kil 



(Pki 

)Bm{k,ki, 



(5) 

■kii; 



In the approach described in [15], Eq. (4) is only a valid 
description of tracer correlations as long as it is indepen- 
dent of the value of the coarse-graining scale Rl. In order 
to confirm this, we need to write Eq. (5) in a form that is 
separable in r and Rl- For realistic bispectra B, this can 
in general only be done approximately. Fortunately, in 
the case of primordial non-Gaussianity where the contri- 
bution in Eq. (5) typically becomes significant for large 
separations r, we can make use of a separation of scales: 
while the Fourier integral over k roughly picks out scales 
of /c ~ l/r, the integral over fci peaks for fci ^ 1/Rl- 
We thus perform an expansion of Eq. (5) in powers of 
fc/fci ~ fci?L. 

We will restrict to local primordial non-Gaussianity 
here for simplicity. As shown in [15], to leading order in 
powers of fc/fci. 



(<5i(l)<5i(2))=4/NLaiC05(r). 



(6) 



Clearly, this contribution is strongly i?L-dependent 
through the factor ct|. The solution introduced in [15] 
is to explicitly account for the dependence of the tracer 
density on the local amplitude of small-scale correlations, 
y* = (<^s/Cs ~l)/2- The corresponding renormalized bias 
parameter 601 then absorbs the term Eq. (6), leading to 
an i?i-independent final expression for the 2-point cor- 
relation. 



The second term here corresponds to the scale-dependent 
bias identified in [10], where 601 quantifies the amplitude 
and is given by the derivative of the tracer abundance 
with respect to a change in the amplitude of primordial 
fiuct nations. 

Unfortunately, Eq. (6) is only an accurate approxima- 
tion for the full expression Eq. (5) on very large scales. 
Taking the expansion in fc/fci to order (fc/fei)^, we will 
show below that 

{Sl(1)SI{2)) = 2/nl [<jU<ps(r) + <Jxx^^s{r)] , (8) 

where (p(x) = — V^(/)(x), and the spectral moment ax,L 
is defined below [Eq. (12) with Eq. (29)]. The crucial 
point is that the second term in Eq. (8) scales differently 
with r and Rl than the first term. Hence, there is no 
hope that it will be absorbed by the single bias parameter 
601 introduced in [15]. In fact, the physical interpretation 
of the two terms in Eq. (8) is quite different: while the 
first term quantifies the uniform rescaling of the local 
small-scale fiuctuations by (j> [Eq. (2)], the second term 
corresponds to a scale-dependent rescaling of small-scale 
fiuctuations by the field </?(x). 



P{k,) ^ [1 + 4/NL^(x)X(fc,)]P(fc,) . 



(9) 



Ur) 



bloUr) 



2bioboi2 f-^i^^^s{r) 



(7) 



The first effect is naturally captured by allowing for a 
dependence of the tracer density on the amplitude of 
small-scale fluctuations. On the other hand, to capture 
the effect of the second term in Eq. (8), we need to al- 
low for a dependence of the tracer density on the shape 
of the power spectrum of small-scale fiuctuations. As 
we will see, such a dependence along with the associated 
renormalized bias parameter is exactly what is needed 
to absorb the i?L-dependence introduced by the second 
term in Eq. (8). 

The purpose of the following sections is to make 
these statements rigorous. Further, we will present all 
derivations for a general bispectrum of primordial non- 
Gaussianity, as the treatment is easily phrased to encom- 
pass the general case. 

The outline of the paper is as follows. In Sec. II, we 
introduce notation and conventions used throughout the 
paper. Sec. Ill describes the squeezed-limit expansion of 
three-point correlators for general (separable) bispectra. 
The scale-dependent bias beyond the squeezed limit is 
derived in Sec. IV, while examples and numerical results 
are presented in Sec. V. We make the connection to 
previous results in Sec. VI, and discuss other sources that 
become relevant in this regime in Sec. VII. We conclude 
in Sec. VIII. 



II. NOTATION 

We assume that the primordial iV-point functions are 
given in terms of the Bardeen potential during matter 
domination, (j). Throughout, we will only deal with the 
statistics of (j) and the initial (linear) density field, scaled 



to some redshift z. Note in particular that the density 
bias parameters correspondingly denote Lagrangian bi- 
ases throughout. The relation between and the linear 
density field at redshift z is written in Fourier space as 



(5(k,z)=X(fc,z)0(k) 



M{k,z) 



2 PT{k)g{z) 
3nraH^{l + z)' 



(10) 



where T(fc) is the matter transfer function normalized to 
unity as fc — > 0, and g{z) is the linear growth rate of the 
gravitational potential normalized to unity during the 
matter dominated epoch. In the following, we will drop 
the argument z in S and Ai since it is not of relevance in 
the derivation. 

Further, we define Myik) — M{k)WYik) where Y 
stands for different filters such as L, s, * which we will 
encounter below. We let P,p{k) denote the power spec- 
trum of (j), and Pm{k) = Ai'^{k)P^{k) the matter power 
spectrum (again dropping the argument z). 

We will need various spectral moments of the density 
field. We define 



af. 



d^k 



Pm{k)W^{k) 



(11) 



for any filter Wy- Further, given a weighting function 
f{k), we define 



2 



d^k 

(2^ 



f{k)P^(k)W^(k). (12) 



It will also be useful to define non-local transformations 
of the density field. Again given a function /(fc), we let 

/li^k 
^^e^^-f{k)Wrik)Sik). (13) 

In particular, this yields 

{Sf^YSy) = crjy ■ 

III. BEYOND THE SQUEEZED LIMIT 



(14) 



We begin with deriving the correlator Eq. (5) and ex- 
panding it in the squeezed limit. More details are pro- 
vided in App. A. Throughout, wc will work to leading 



order in the dimensionless amplitude of non-Gaussianity 
/nl- At this order and for the models we consider, the 
only relevant A^-point functions arc the power spectrum 
and bispectrum. 

Expressed in terms of the bispectrum of primordial 
perturbations B^, defined through 

(</'k0ki0k,) = {27Tf6Dik + ki + k2)B^(k, ki, k2) , (15) 

the term wc arc interested in becomes 



{Sl{1)SI{2)) = 



d^k 



^''■''MLik) 



d^ki 



(27r)3 "^^ ' J (27r)3 

xXi(|k + ki|)S^(fc,A:i,|k + ki 



ML{ki) 



(16) 

Our goal is to obtain an expression that is separable in 
r and R^. Neglecting the i?L-dependence introduced 
through 7Wi(fc), which is irrelevant if r ^ Rl (see 
Sec. VII), separability in r and R^ is equivalent to hav- 
ing an integrand separable in k and ki. As a necessary 
prerequisite, we assume that the bispectrum B^ is given 
in separable form, 

B0(fcl, fc2, fcg) - 5] A„ [FW(fci)F^2)(^^)^^3)(fc^) 



-|- 5 perm. 



(17) 



where Aa are constants and the 6 permutations guar- 
antee the symmetry of B^ in its arguments. The sum 
a runs over however many terms are necessary to ade- 
quately approximate the bispectrum in separable form. 
We further define 



F^{k)=ML{k)F^(k). 
Eq. (16) then becomes 



{k)e 



ik-r 



(18) 



(19) 



j3 L 



We now expand the ki integrand in powers of k/ki up to 
second order. As shown in App. A, 



{S,{l)Sl{2)) =Y^A^y ^^Fi'\k)e^^- 

d^ki 



(27r)3 
2kl 



2kl 



2Fi'\k,)Fi^\k,) + ^-^Fi'Hk,) (1 - M^)fci/^(^)(fci) + M^fc?F:(3'(fci) 



Fi^Hk,) (1 - f^')k,F^^'\k,) + ^?klF'^^\k{, 



+ {(123) ^ (231)} + {(123) ^ (312)} , 



where primes denote derivatives with respect to fci . This 
expression is now in the desired separable form. It is vahd 
up to terms of order fc^/fcf as the cubic terms drop out 
just hke the hnear terms did. We can make this result 
more obvious and compact by introducing the notation 



(2^3 



^SW^ / 7!^>'L(fc)F«(fc)e*- 



i^lL(r)^ J ^s-^L(k)eFiHk)e^'-'' 



'aL 



1 

2^ 



dkiklMl{k,)F^\k,)Fi^Hk,] 



(20) 
(21) 

(22) 



'XaL 



127r2 



dfci{F«(fci) 

X "2fcii^(^")(fci) + fc2/'^'W)(fci) 
+ W^ (.?)}• (23) 



Note that for all spectral moments, a^^"^^^ — cr^'^*^ and 
that the derivatives in Eq. (23) act on both Fa'{ki) and 
ML{ki)- With this, we obtain 



{5L{i)5im) ^Y.M sciW'^'f^ +4'IlM'^ 



2(23) 
^XaL 



+ {(123) ^ (231)} + {(123) ^ (312)} 

(24) 

where the second line denotes the two remaining cyclic 
permutations. Apart from the residual i?i-dependence 
in i,aL, S.v^aLj this expression is fully separable in r and 
Rl as desired. The expansion in k/ki can of course also 
be taken to higher order if necessary. We will discuss the 
necessity of this below. 

Given the appearance of ct^^ and dxaL i^ Eq. (27) , the 
non-Gaussian correction to ^/j is strongly i^^-dependent. 
In the next section, we will address this issue through 
renormalized bias parameters. 



A. Primordial non-Gaussianity of the local type 

Eq. (24) is general, but somewhat abstract. In order 
to clarify its physical significance, we now specialize to 
the case of local primordial non-Gaussianity, for which 
the bispectrum is given by 

B^(ki,k2,k3) -2/NL[P0(fci)/'0(A:2) + 2pcrm.]. (25) 

Since B^ is already separable, wc only have one term in 
the sum Eq. (17), with Aa = A = /nl, and 

F'^^Hk) = F''^\k) = P4k); F(3)(fc) = l. (26) 

The permutation CaLi''")'^aL ^^ ^Iq. (24) is suppressed 
(for scale-invariant P^) by {k/kiY relative to the other 



two identical permutations. In keeping with our treat- 
ment up to (fc/fci)'*, we will retain the leading contribu- 
tion from this term. We then obtain 

+ 2h^.T'k^:i{r), (27) 

where we have defined 

^<^sAr) = j^ML{k)P4k)e^^- = {HI)Sl{2)} (28) 
^<fSr.{r) = ^MUk)k'P4k)e^'^-'^ = Ml)SLi2)) 
Xik) = ^{M-i(fc) [2kM'Ak) + eM'[{k)] (29) 
+ (P^Ml)^' [2fc(P0ML)'fe + k\P^ML)t] } , 

and (t|, (Jxl ^^'^ defined through Eqs. (11) (12). Here 
we have introduced (/^(x) = — V2(/)(x). The contribution 
from the third permutation is given by 



d'i«-/(g3-^^(fc> 



2(12) ^ 
'aL 



27r2 



dkiklP4,{ki)Pm{ki)Wl{ki). (30) 



The first (leading) term in Eq. (27) is well known and 
agrees with that derived in [15]. This term can effectively 
be described as a rescaling of the local density field. 



,5(x) ^ [1 + 2/nl'^(x)] (5(x) , 



(31) 



which leads to Eq. (2). On the other hand, the second 
(subleading) term in Eq. (27) can be seen as coupling the 
Laplacian of (f) to the density field. However, this cou- 
pling is scale-dependent, i.e. when performing a Fourier 
transform in a local patch where ip can be considered 
constant, we have 



<5(k) ^ [1 + /NLA(fc)(^(x)] S{k) . 



(32) 



This effect is of course suppressed with respect to cj) in 
the large-scale limit. In Fourier space on large scales 
M.{k)P^{k) ex k~^Pjn{k), so that ^.psir) grows with re- 
spect to ^L(r) on large scales, while £,ips.L{''') ^ Ci(^) (the 
nontrivial transfer function leads to departures on scales 
smaller than about 100 h^^ Mpc). 



IV. SCALE-DEPENDENT BIAS 

As shown in [15], one can absorb the leading term oc af^ 
in Eq. (27) by introducing a dependence of the tracer 
density on the local variance of the small-scale density 
field. We define the small-scale density field as the local 



fluctuations around the coarse-grained fleld Sl' 

,5,(x)=(5,(x)-<5i(x) (33) 

= J d^y[W,{^ -y)-WL{x- y)]5(y) 



^3u 

W^,(fc)J(k)e'''•^ 



(27r)3 



(34) 



where we have introduced a fixed small smoothing scale 
-R*. Thus, we write 



nft[5L(x)] ^n,i[(5L(x),y*(x)] 



(35) 
(36) 



where the subscript * refers to the smoothing scale i?*, 
(?/*) = 0, and the factor of 1/2 is included to obtain 
expressions which conform to standard convention. A 
very similar derivation to that of Eq. (24) yields (App. A) 



2(23) 
as 
2 



+ 2 perm. I , (37) 



where a subscript s indicates that the filter function Wl 
should be replaced with Ws in Eqs. (22)-(23). In the case 
of local primordial non-Gaussianity, this again reduces to 



1 



(<5L(l)y.(2)) - 2/nl Us.ir) + T^^^sjr) 



1 2(12) 

I- ^{3) 0'as 
2?aL 9 



'Xs 



(38) 



While the renormalization of the -R^-dependence in 
Eq. (24) works for the leading term as will see now, we al- 
ready notice two issues with Eq. (38) in comparison with 
Eq. (27). First, the second and third term in Eq. (38) 
depend on R^, since Wg is defined with respect to Wl 
[Eq. (34)]. Thus, the final expression cannot simply in- 
volve ((5L(l)y*(2)). Second, there is no reason why the 
relative magnitude of the three terms in Eq. (38) should 
be the same as for the three terms in Eq. (27), so that 
we cannot expect all terms in Eq. (27) to be absorbed by 
((5l(1)j/*(2)). These issues will be resolved below. 



A. Bivariate bias expansion 

We begin by re-examining the bivariate bias expansion 
of [15] (which in the local case at lowest order is equiva- 
lent to what was presented in previous papers [18, 19]). 
Throughout this section and the next, we will restrict to 
a single separable contribution a and a single permuta- 
tion (123) of Eq. (24) [note that this corresponds to two 



permutations in Eq. (17)]. The complete prescription for 
tracer clustering, presented in Sec. IV C below, will then 
involve a sum over the contributions from different a and 
permutations. 

In terms of the "bare" bias parameters c„m, that is, 
the coefficients of the Taylor series of rih in 6^, y^, the 
tree-level expression for the tracer correlation function in 
the bivariate PBS expansion is [see Eq. (109) in [15]] 



C^°M - -^\ cUL(r)+cioC2o{5L{l)Sl(2)) 



1 



1 



■ 2cioCoiAq 



2cioCoi(<5L(l)y,(2)) 



(39) 



t(i) 



2(23) . f(l) 



2(23) 



'^'iaLVJ^aL +^W^aLyl"XaL 



2(23) -, 2(23) n 

(1)^.\'^"'' _L^c(l) (^\'^Xas 






where 



A^- E 



^nm /cn^.ni 



1 'm ' 



{Sly: 



(40) 



n,m— 



Let us focus on the leading terms in the squeezed limit 
first. We need to generalize the definition of the rcnor- 
malized bias parameters bj\jM to the case of the general 
separable bispectrum Eq. (17). 

Consider a rescaling of the Fourier-space density field 
given by 



S{ki)^5{ki)[l + eF{ki 



Fiki) 



Fc 



(2) 



{k,)F^^\k,) 
P^{ki) 



(41) 



We then define the bivariate bias parameters Iinm as the 
response of the mean tracer density to a change in the 
background density (at fixed proper time) and the rescal- 
ing of the initial density field under Eq. (41): 



iiNM — 



1 



a^+^^(n„)D,, 



{nh)D=o,e=i) dD^de 



N P)cM 



(42) 



D=0,e=0 



«(23) 



Note that b^M should really be denoted &jvV here, since 
the definition refers to the specific rescaling in Eq. (41). 
However, since we are only dealing with a single a and 
permutation here and in Sec. IV B, we will drop this des- 
ignation for notational clarity. Note further that the di- 
mension of 6oi is such that the contribution to the corre- 
lation function is dimensionless. In particular. 



dim(&oi) =dim[F(fci)] = 



Mpc^ 



dim(Fi^^) 



(43) 



where we have assumed that the Aa are dimensionless 
so that dim(F^^Vi^Vi^^) = Mpc*^ [Eq. (17)]. Under the 
rescahng Eq. (41), we have [using the definition Eq. (13)] 



(5l(x) -^ (5l(x) + e (5f,l(x) 



2/*(x) -^ y*(x) + ■^ds(x.)dF.s{x) + TT^Sl (x.) 



(44) 



Inserting the expression for 6oi [Eq. (46)] yields 
1 



a™°™(^) = T72 ^lo^i(^) + 2ciocoi^^A„Ci:(^) 



2(23) 
O'as , ,(1), 



M- 



2(23) . ,(1) 



+2cioc2o<rA„ei:M . (49) 



Note that 



((5y(x)(5F,y(x)) = cr|_y = a]!^y^ 
for y = L, s, . . . . Using Eq. (35), we obtain 



If the renormalized bias 601 defined with respect to 
the transformation Eq. (41) properly removes aU Rj^- 
(45) dependence of the tracer 2-point function, Eq. (48) 
should match the bare bias expansion Eq. (39). Let us 
thus take the difference: 



^nh 



)iD,e)^{n,{0))J2- 



?, ' T77, ' 



{Sl+sSf,l)Sl + D 



;SsSf,s 






To the order we are interested in, we only need boi, which 
is given by 



"01 



\l ^ — ^ r)}m} 



N 



X U{5F,L5l-'yT) + "^^{Sl^sSp^syT-^)] 



2(23) 



NV''- a^. 



2(23) 



C20Cr„i 



0{5\f^^) 



(46) 



Following the same reasoning as in [15], our guess for the 
tree-level tracer correlation function written in terms of 
renormalized bias parameters is 



er°™(r)=6?oaM + 2foio&oi 



2(23) 



{5L{l)y.{2)) 



(47) 



2(23) 

+ 2feio&oiA„Ui^i(r) + ^e^\)„^(r)''^- 



2(23) 



where the prefactor in front of ((5L(l)y*(2)) takes into 
account that the tree-level relation between 5oi and cqi 
is 601 = coiCas /cg. In the second line we have used 
Eq. (37). The fact that renormalization is not success- 
ful at subleading order can already be seen from the last 
term in this equation. This depends on Rl through the 
spectral moments, since the kernel Ws depends on R^. 
Thus, we assume that a proper renormalization will ab- 
sorb this term as well, and write our updated guess as 



^--™(r) = bUdr) + 2610601 Aoeii!(0 . 



(48) 



iH{r)^i 



rcnorm 
h 



(r) 



(50) 



^aCv2„i(?-)cio 



2(23) 
C20<yxcL 



coi 



2(23) 
^Xas 



All terms here are sub-leading in the squeezed limit, that 
is, at leading order in the squeezed limit the renormaliza- 
tion works in the same way as shown in [15]. However, 
when going beyond the squeezed limit for general primor- 
dial non-Gaussianity, the description of the tracer density 
as a bivariate function n/JfJ^, y»] is not sufficient. 



B. Trivariate bias expansion 

We are seeking a third local parameter (besides 5l and 
y■^,) which the tracer density in general depends on, and 
which will allow us to absorb the i?j,-dcpendcnt terms in 
Eq. (50). In the local model, the residual i?L-dependence 
that we encounter when going beyond the squeezed limit 
is induced by the fact that in this case, the effect of long- 
wavelength modes on small-scale modes is not simply to 
rescale them uniformly, but to also change the shape of 
the small-scale power spectrum [Eq. (29)]. More gener- 
ally, small-scale modes are rescaled differently as function 
of the correlation scale r, as the relative importance of 

the terms given by Cali''") ^^^ w^ail'') changes. Thus, 
we have to allow for a dependence of the tracer density 
not only on the local variance of small-scale fluctuations, 
but also on their shape. As a proxy for the power spec- 
trum shape, we will use 
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d\nR, 



-y*(x). 



(51) 



This choice is useful because transformations of the den- 
sity field of the form 



5l + D- (5(x) ^ [1 + t] (5(x) 



(52) 



leave /x*(x) invariant. Using Eq. (51) and Eq. (37) we 
obtain, again for one of the three cyclic permutation of 



the contribution a, 



{SL{l)^^.m) 




(53) 



For local non-Gaussianity, (7^^ = a'^, so that the first 
term drops out, and we obtain (restoring the number of 
permutations) 
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Let us now adopt /i* (x) as third local parameter: 

rih [Sl (x) , y* (x)] -^ uh [5l (x) , y* (x) , ^* (x)] . (54) 

We again consider the scale-dependent transformation 
Eq. (41), and in addition a further, independent scale- 



dependent transformation of the density field 

6{\i) ^ [1 + Lf{k)] ,5(k) . (55) 

Here we leave the function /(fc) free for the moment. This 
leads to 
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In keeping with the leading order treatment in /nl, we 
will only consider the first derivative of the mean tracer 
density with respect to t. Thus, we can drop the i^ terms 
in Eq. (56). We define 
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where the mean tracer number density is given by 
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We thus obtain 
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(60) 




where we have expanded to order 5\. Note the non- 
trivial coefficient multiplying cqoi in Eq. (60). We thus 
have to divide feooi by this coefficient when writing down 
the contributions to £^h that correlate /i* with other per- 
turbations. 

We now update our ansatz for the renormalized tracer 
correlation [Eq. (48)] to the trivariate case. Since we 
found in the previous section that the bivariate case 
(neglecting the dependence of Uh on /x*) successfully 
produced an i?^ -independent expression for the leading 
squeezed-limit contribution, we expect that only the sub- 
leading contribution of ((5^(1)//* (2)) [Eq. (53)] will con- 
tribute. That is, we expect the leading contribution in 
Eq. (53) to be absorbed by 6oio- Thus, our expectation 



C'\r) = bj,,^L{r) + 2h,obomA^Cii{r) 
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The last term is in general only guaranteed to be Rl- 



independent if (JxaJ = a'i g, that is, if 



klPm(h)f{k,) = lFi'\h) [2k,F:J^'Hk,) + klF'^^'Hk,) 



+ (2) ^ (3) . 



(62) 



Thus, we will fix this as our choice of f{k), and in the 

following derivation set cTxaY ~ '^I'Y^ where Y = s,L. 
Note that the dimension of 6ooi is equal to the dimension 
of / and is given by 
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that is, the dimension of 6oi times Mpc^. The bare ex- 
pansion then becomes [using ajy — '^f y from Eq. (45)] 
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On the other hand, inserting Eqs. (59)-(60) into Eq. (61) yields 
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We thus find that the trivariatc renormalized expression, 

1,2 
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(64) 
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(66) 



has successfully absorbed all i?L-dependence from the bare expansion including subleading terms in the squeezed limit 
(that is, apart from the residual dependence through ^l, ^l^£, Cv^aLi which is negligible on large scales). 



C. Total contribution 



We are now ready to derive the full expression for the two-point tracer correlation including subleading corrections in 
the squeezed limit for the general separable bispectrum Eq. (17). For each a, and each of the three cychc permutations, 
we define general trivariate bias parameters through 
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where e("'*J) , i(".*i) are scale-dependent rescalings of the density field given by 
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F«(fci) \2k,F^^Hk,) + klF^^^'Hh)] + (z) O (j) 



(68) 



Note that b^j^l ~ ^nml- Then, the renormalizcd tracer 2-point correlation at tree level is given by 



d"M=6?ooeL(r) + 26iooE^" 



iT^^^^kr) + ^4f ^41l(^) + {(123) ^ (231)} + {(123) ^ (312)} 



(69) 



Note that when going beyond tree level, one in general also has to include mixed contributions simultaneously 
involving rescalings Eq. (68) for different a and/or permutations. Eq. (69) is accurate up to terms of order (fc/fci) in 
the squeezed limit. We will make this statement more precise in the following sections. We will discuss the expected 
relative amplitude of feoio a-nd feooi below in Sec. VB. Note again that both 6oio ^nd 6ooi are in general dimensionful 
[Eq. (43), Eq. (63)]. ^ 



D. Tracer correlations in Fourier space 

In Fourier space, the tree-level result Eq. (69) can be 
phrased in terms of a generalized scale-dependent bias: 
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Thus for each contribution a in the separable bispectrum 
Eq. (17), and for each of the three cyclic permutations, 
we have in general two independent contributions to the 
scale-dependent bias. The leading term in the large-scale 
limit {k — >■ 0) scales as S{k), while the subleading term 
scales as k'^S{k). 

Eq. (70) now allows us to estimate the relative magni- 
tude of the subleading term. For each a and permutation. 
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On the other hand, the contributions from different a 
and permutations (ijk), (Imn) scale as 



I3ir. 
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(73) 



and similarly for the subleading term. Since Eq. (69) 
and Eq. (70) are derived neglecting terms that, for each 
contribution a and permutation, are suppressed by k* 
relative to the leading contribution, only contributions 



which are suppressed by less than k^ in Eq. (73) rela- 
tive to the overall leading contribution [i, a) should be 
included in Eq. (70). 

Note further that we have not included the stochas- 
ticity in Eq. (70) which in some models of primordial 
non-Gaussianity can become important on large scales 
as well [20, 21]. 



EXAMPLES AND NUMERICAL ESTIMATES 



A. Local non-Gaussianity 



As before, for local non-Gaussianity we consider the 
two leading cylic permutations as well as the leading con- 
tribution from the third permutation. Thus, 
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Here, the bias parameters are derived with respect to the 
rescalings 



.loc 

"oio 

tIoc 

"ooi 

,(12)loc 

"oio 



,5(ki)^ [l + e]<^(ki) 
<5(ki)^ [l + eX{k,)]S{k,) 
<5(ki)^ [l + eP4fci)]5(ki) 



(75) 



10 



where X{ki) is defined in Eq. (29). Note that the leading 
term agrees with [15]. In Fourier space, this becomes 
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FoUowing Sec. IV D, we can estimate the relative mag- 
nitude of the subleading correction and the contribution 
from the third permutation relative to the leading term 



^"subleading (''') 

A5(12)loc(fc) 
A&!°^ding(fc) 



Lloc 

"001 7 2 
oUoc 
^"010 

,(12)loc 

"oio 



5[,To 2F^(/c) 2Asb'i 



,-3i,(12)loc 

'^0 "OlO jl3 

oc 



010 



(77) 



where in the last line we have assumed a scale-invariant 
spectrum of (j) with amplitude ,4s defined at the pivot 
scale fco, 
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Thus, the subleading term is suppressed by a factor of fc^ 
with respect to the leading term, while the contribution 
from the third permutation is suppressed by k^. 



B. Universal mass function 

In order to quantitatively assess the importance of the 
subleading term in Eq. (69) and Eq. (70), we need to 
estimate the magnitude of 6ooi- For dark matter halos 
this can be done accurately through N-body simulations 
with modified initial conditions following Eq. (68) [or 
Eq. (75) for local non-Gaussianity]. However, a detailed 
comparison with N-body simulations is beyond the scope 
of this paper. Instead, we make use of a generalization 
of the universal mass function prescription as discussed 
in Sec. IV E of [15]. We write the mean abundance of 
tracers as 



FIG. 1: Contributions to the scale-dependent bias from local 
non-Gaussianity (/nl = 1) for halos with M = 2-10" h~^ Mq 
at 2 = (6ioo — 0.44) and assuming a universal mass function 
from the Sheth-Tormen prescription [22], scaled by (k/Ho)^ 
to yield a scale-independent value on large scales. The dashed 
line shows the leading term, the solid the leading plus sublead- 
ing (order (fc/fc*)^) term, while the dash-dotted line includes 
the order (fc/fe«)^ term. The red long-dashed line shows the 
prediction from the conditional PS mass function (Sec. VI). 



smoothed on a single scale i?*, as well as its derivative 
with respect to scale. 

Let us again consider an individual contribution a and 
permutation (123). Under the rescaling Eq. (68), cr* 
transforms to lowest order as 



nh =nh {p,cr^,, J*) 
_ diner* 



d\nR^ 



(79) 
(80) 



where ct* is the variance of the linear matter density field 
on scale i?*, and i?* is related to the mass M* through 
M* = 47r/3pi?2- The Jacobian J* is present in order to 
convert from an interval in ct* to a mass interval. In this 
approximation, n^ is given as a function of the mean den- 
sity of the Universe and the variance of the density field 
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where f{k) is defined through Eq. (62). As shown in [15] 
(see also [12]), the scale-dependent biases are then given 
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Here, Uqiq is the leading scale-dependent bias param- 
eter for local primordial non-Gaussianity, for a tracer 
following Eq. (79). We have assumed that the tracer 
density scales linearly with the Jacobian as expected 
physically. For such tracers, the leading and sublead- 
ing bias parameters quantifying the response to gen- 
eral non-local non-Gaussianity are thus directly related 
to the leading bias parameter for local non-Gaussianity. 
If we further specialize to a universal mass function, 
fih = fihip,!^ = 5c/cr*,J*), then 6o°o == &ioo(5c (recall 
that &100 is the Lagrangian bias). 

The precise magnitude of the bias coefficients 6oio, ^ooi 
depends on the exact rescaling Eq. (68). Generally, if the 



main contribution to CTc, 
fc ~ fc*, then 



'/,* 
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comes from Fourier modes 
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Note that depending on the shape of the rescaling fc» 
does not necessarily have to be of order 1/i?*. Thus, for 
tracers following Eq. (79), Eq. (72) simplifies to 
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More generally, Eq. (79) implies that for scale-free bis- 
pectra for which Sa (k) (x fc"° , fc* is the only scale in- 
volved, so that other contributions as in Eq. (73) will all 
be suppressed by powers of k/k^,. Conversely, k ^ k^ in- 
dicates the breakdown of the perturbative expansion in 
the squeezed limit. 



C. Numerical results 

Fig. 1 shows the leading and sub-leading scale- 
dependent bias contributions for local non-Gaussianity 
assuming tracers following a universal mass function. 
Specifically, we assume a halo mass Af = 2 • 10^^ h~^ Mq, 
so that 6g°'5 — biooSc — 0.44. We have multiplied the 
scale-dependent bias by {k/Ho)"^ in order to obtain a 
weakly scale-dependent result. Note that even the lead- 
ing term has a residual scale dependence due to the trans- 
fer function contained in A4{k). We also show the con- 
tribution A&(^^''°'^(fc) from the third permutation of the 
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FIG. 2: Leading and sub-leading contributions to the scale- 
dependent bias from local non-Gaussianity, as in Fig. 1, but 
for different masses at z = 0. The curves shown correspond 
to, from top to bottom, M = 2 ■ 10^'', lO", 5 • lO", 2 ■ 10^^ 
10", 5 • 10l^ 2 ■ 10^^ and 10^^ h'^ Mq, respectively. 



local bispectrum [Eq. (77)]. The higher order contribu- 
tions become important as fc > 0.05ft.Mpc~^. Note in 
particular that the term Afe'^^^^'°'^ grows rapidly towards 
smaller scales. We will discuss this issue in the context 
of the relation to previous approaches in Sec. VI. 

Fig. 2 shows the leading and leading-|-subleading con- 
tributions to the scale-dependent bias for a range of 
masses from 10^^ to 2 • 10^'' h~^ Mq. Clearly, the typical 
scale at which the subleading correction becomes impor- 
tant does not depend sensitively on the mass. Gener- 
ally, the correction is more important at lower masses, 
specifically around M^ where the Lagrangian bias 5ioo 
vanishes. For vanishing feioo, the leading term vanishes, 
whereas the subleading term does not disappear entirely 
due to the non- vanishing derivative of In a^ with respect 
to In a'i . Since aj^ * scales more weakly with R^ than a^ , 
this derivative is less than one leading to a suppression of 
the scale-dependent bias when including the subleading 
correction. 

In principle, a general tracer could lead to very differ- 
ent numerical results, i.e. much larger or smaller sublead- 
ing corrections. However, given the accuracy of universal 
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mass functions of 10 — 20% at least for dark matter halos, 
we expect the magnitude of the corrections as shown in 
Fig. 2 to be typical. 



VI. RELATION TO PREVIOUS RESULTS 

We now make the connection to previous results on the 
scale-dependent bias for general non-Gaussianity beyond 
the squeezed limit. In Desjacques et al. [12], the scale- 
dependent bias was derived by applying a conditional 
mass function approach to the Press-Schechter (PS) mass 
function [23]. The non-Gaussianity was taken into ac- 
count by applying an Edgeworth expansion to the Gaus- 
sian PDF of the density field. As shown in [12], the 
scale-dependent bias defined through Eq. (70) is in this 
case given by 
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Here, a subscript * denotes filtering with a tophat of 
radius R^, the Lagrangian radius corresponding to the 



mass scale of the tracer as in Sec. VB (this was denoted 
as Rs in [12]). Note the derivative with respect to IncrJ 
rather than Incr* as written in [12]. Here, the function 

Ti (fc) is given by an integral over the bispectrum, 

X B0(k,ki,-k-ki), 
(87) 

where M.^{k) = M{k)W■^,k). The function J"i ' can then 
be directly related to the squeezed three-point function 
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again up to order {k/k^)"^, see for example the derivation 
leading up to Eq. (A7). In Fourier space, this relation 
becomes 
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Eq. (86) then becomes 



^'^^^'^-pMMm^^A^'^'^^ 



E^-i'^^'^w 



hi5c 



biSc 



d 



dlnal 



diner? 



2(23) -, 



2(23) , 



biSc 



biSc 



din 0-2 



diner? 



2(23) 



2(23) 
cr? 



2 perm. > , 



2 perm. 



(90) 



r 



using the definitions after Eq. (70). Comparing with 
Eq. (70), we can now read off the bias parameters 



,(a,23) 
^10 



, (a,23) _ 

^ol ~ 



6i(5, + 2-^ 
'b5 nl"^^"" 


2(23 


\ 
-I 
J . 




2(23) 


2(23) 




<T? 


biS, + 2 -—- 
\ dlr 


2(23) 
Xa* 


y 
-1 


2(23) 
'^Xa* 



(91) 



(92) 



We see that both the leading and subleading bias pa- 
rameters derived from the conditional PS mass function 
agree with those expected from a general universal mass 



function [Eqs. (82)-(83), with &q°q — biooSc]- Fundamen- 
tally, this is a consequence of the fact that in the Press- 
Schechter approach, as in general for universal mass func- 
tions, there is only one scale _R* that enters the descrip- 
tion of tracer statistics. We thus expect this result to 
hold at higher order in fc/fc* as well. 

Since Eq. (86) does not involve a perturbative expan- 
sion in the ratio of wavenumbers k/k^ , it can serve as 
a useful guide as to where this expansion breaks down. 
Fig. 3 shows the residuals when including the leading and 
subleading terms from Eq. (76). We see a residual which 
scales as k^ for small k. When including the term from 
the last line of Eq. (76), we see that the residuals become 
even smaller and scale as k as expected. Some numerical 
artefacts are visible when the residuals become of order 
10~^ or smaller. These are due to the numerical deriva- 
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tive performed when evaluating Eq. (86). We can in fact 
perform a rough estimate of the expected correction or- 
der [k/k^.Y in Eq. (76), via 
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Note that due to the logarithmic divergence we need to 
introduce a low-fc cutoff in cr^4*. This is likely to be 
an artefact of the universal mass function prescription, 
where the leading effect of a change in the small-scale 
power spectrum shape on the tracer density is given by 
this formally divergent spectral moment. In reality, trac- 
ers will have a finite response to such a change. Here 
we choose kmin = 0.01/iMpc~^, corresponding to the 
turnover in Pm{k). Eq. (93) is in any case only to be seen 
as a very rough estimate. This contribution is shown as 
dotted line in Fig. 3, making clear that the residual, after 
taking into account all terms in Eq. (76), indeed scales 
as fc'* on large scales. 

Fig. 3 shows that for the local model, the order {k/k^,)^ 
correction becomes comparable to the lower order correc- 
tions at fc ^ 0.02 h Mpc^ , signaling a breakdown of the 
perturbative expansion there, even though the fractional 
deviation from the full result Eq. (86) when including 
terms up to order (fc/fc*)'^ remain at 10% or less all the 
way to fc-- 0.1/iMpc~^ 

In summary, the conditional PS mass function results 
derived in [12] are consistent for tracers following a uni- 
versal mass function, in the sense that they match the 
results from the general renormalization approach when 
restricted to universal mass functions. Note that this 
holds once the polynomials in Sc/cr^ are replaced with 
bias parameters, as described in [12]. In this context, the 
key advantage of the PS approach is that it sums over 
all powers of fc/fc* , without relying on a perturbative ex- 
pansion in this parameter. 

On the other hand, realistic tracers will not simply de- 
pend on the variance of the density field on a single scale, 
thus breaking the relation between the bias parameters 
^001 J ^010 1 and 6ooi. Furthermore, there are other scale- 
dependent biases which contribute at the same order as 
the subleading correction 5ooi. We will turn to this issue 
next. 



VII. SCALE-DEPENDENT BIAS BEYOND THE 
LARGE-SCALE LIMIT 

We have seen that beyond the squeezed limit, there is 
a subleading correction to the scale-dependent bias from 
primordial non-Gaussianity that scales as k'^ relative to 
the leading term. In addition to this correction however, 
we expect two additional contributions that are leading 
order in /nl (i-e., in the primordial bispectrum), and 
have the same scaling with k. 
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FIG. 3: Fractional difference between the contributions in 
Eq. (76), evaluated for a universal mass function, and the re- 
sult for the conditional PS mass function Eq. (86) from [12] 
for local primordial non-Gaussianity. The black line solid line 
shows the residuals when including the leading and sublead- 
ing (order (fc/fc*)^) contributions, while the green dash-dotted 
line also includes the order (fc/fc*)'^ contribution from the last 
line of Eq. (76). The dotted line shows a rough estimate of 
the order (fc/fc,)* correction (see text). We have again as- 
sumed M = 2 • lO" h~^ Mq and z = 0, although the results 
are essentially independent of the mass. 



First, as shown in [15, 24], non-locality in the forma- 
tion of tracers generically induces a dependence on the 
curvature of the density field, leading to a contribution 
of 



Ur) ^ b^2sboio{^HL{l)y*m 



(94) 



This contribution will also serve to absorb the residual 
i?L-dependence present in £,^[{r), ^^i^i^ir) in a similar 
way as discussed in [15]. In Fourier space, Eq. (94) cor- 
responds to a contribution to the scale-dependent bias of 
the form 



Phjk) 

P^[k) 



D2&v25fc2A6ioad(fc), 



(95) 



where Afoioad is the leading contribution to the non- 
Gaussian scale-dependent bias from Eq. (71). If Lg is 
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the scale of non-locality of the tracer (in terms of its de- 
pendence on the matter density), then b^2s ^ i^, so that 
this additional contribution scales as (kLg)'^ relative to 
the leading term. 

Throughout the discussion of the non-Gaussian case 
in [15] and here, we have assumed that the tracer den- 
sity is a purely local function of y.^ , the parameter which 
quantifies the amplitude of small-scale fluctuations (in 
this paper, we have also introduced a local dependence 
on /i* — dy^t/dliiR^,). In general, however, one also ex- 
pects that tracers depend on the amplitude of small-scale 
fluctuations in some finite region of size Ly. Then, in 
straightforward analogy with the density case, we also 
need to allow for a bias h^2y with respect to V^y*, where 
h-^2y ~ L^. That is, strictly speaking we need to gen- 
eralize each of the rescalings in Eq. (68) to be spatially 
dependent, e -^ ex^. Schematically, this leads to a con- 
tribution to the tracer correlation of 

Ur) D b,oobv-y{SL{l)V\,m , (96) 

which, in Fourier space, becomes 

Phjk) ^ ,, , ,2 A6lead(fc) 

— -TTT D 2bioabv2yk — . (97) 

This contribution thus scales as (kLy)'^ relative to the 
leading term. 

Typically, one might expect Ls ^ Ly r^ R^,, where i?* 
is the Lagrangian radius of the region that collapses to 
form the tracer. On the other hand, the length scale 
we found for the subleading terms in the squeezed limit 
is 1/A:* ~ 50/i~^Mpc, suggesting that this correction is 
somewhat more important than the other contributions 
described in this section for typical tracers for which 
i?* = 1 — 10/i~^Mpc. However, the value of fc* de- 
pends on the specific type of non-Gaussianity considered, 
and in general all three subleading contributions to the 
scale-dependent halo bias can be comparable in magni- 
tude. Thus, if one of them is included (even implicitly as 
for example in the conditional PS mass function result), 
then all of them should be included for consistency, unless 
one can show that the subleading contribution dominates 
over the other contributions described in this section. 



VIII. CONCLUSIONS 

We have derived the subleading contributions to the 
scale-dependent bias Ab{k) of large-scale structure trac- 
ers for a general separable primordial bispectrum. The 
leading contribution is given by the scaling of the bis- 
pectrum in the squeezed limit, linife^o B(t>{k, kg, jk,, + k|), 
and the scale-dependence of the subleading contribution 
is suppressed by a factor of k"^ relative to this term. 
This subleading contribution is important to quantify, 
since it tells us at which k the usual squeezed-limit re- 
sult ceases to be accurate. For local non-Gaussianity 
and tracers following a universal mass function, we found 



that this happens at fc ^ 0.02 h Mpc""'^, although the first 
two leading terms provide an excellent approximation to 
the result from the conditional PS mass function up to 
k < 0.1/iMpc~ . Our approach is independent of any 
assumptions on the tracers apart from a finite scale of 
non-locality. 

Throughout, our results have been phrased in terms of 
the three-point function of the primordial perturbations, 
which is the standard result of computations performed 
for particular inflationary models, thus allowing for a di- 
rect application of the results of this paper to models of 
inflation. In contrast, several previous papers [11, 14, 19] 
have employed a fictitious Gaussian field mapped to the 
physical field via a quadratic kernel. In this latter ap- 
proach, the effect on large-scale structure tracers is medi- 
ated by modulated spectral moments of the density field, 
e.g. <j1\tf,. This approach is complicated by the fact that 
the kernel is not uniquely determined by the bispectrum, 
so that additional constraints need to be imposed [14]. 
However, the derivation in this paper can also be applied 
to the kernel approach in a straightforward way. Note 
that while the bispectrum specifies the kernel uniquely in 
the squeezed limit [11] (for non-divergent kernels), this 
is no longer the case when including subleading terms. 
Thus, different kernels which yield the same bispectrum 
are expected to make different predictions for the sub- 
leading contribution to the scale-dependent bias. We 
leave this issue for future work. 

In addition to the subleading contribution in the 
squeezed limit, we have pointed out two other contri- 
butions that will contribute at linear order in /nl and 
with the same tree-level fc-dependence (Sec. VII). These 
contributions should be included (or at least carefully 
considered) when putting constraints on non-Gaussianity 
using large-scale structure statistics on intermediate and 
small scales, i.e. for k > 0.05 h Mpc~ . 

It is straightforward to extend the squeezed limit ex- 
pansion presented here to higher order in k. In that case, 
one needs to add a dependence of the tracer density on 
another property of the density field further quantifying 
the sensitivity to the amplitude of small-scale fluctua- 
tions as a function of scale. One possible choice would 
hed^y,/d{\nR,f. 

Another straightforward extension is the inclusion of 
higher primordial A''-point functions. For example, in 
the presence of a primordial four-point function both the 
linear and quadratic bias become scale-dependent [12, 
25], and one needs to take into account the dependence 
of the tracer density on the skewness of the density field. 
No conceptually new issue arises, and the calculation will 
closely follow the one presented here. 

Finally, we have shown how the coefficient of both lead- 
ing and subleading terms in the scale-dependent bias for 
a general bispectrum can be derived for tracers identi- 
fied in N-body simulations, by running simulations with 
modified initial conditions [Eq. (68)]. This will allow for 
a precise test of the accuracy of the universal mass func- 
tion prediction for dark matter halos in the context of 
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non-Gaussian halo bias. We leave this for future work. 
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Appendix A: Derivation of Eq. (24) and Eq. (37) 

In this appendix we derive the next-to- leading order squeezed- limit expressions Eq. (24) and Eq. (37). As described 
in Sec. Ill, wc assume that the bispectruni B^ is given in separable form. 



B^{ki,k2,k3) =^ 



+ Fi'\k,)Fi'Hk2)Fi'Hh)+Fi'Hki)Fi'Hh)Fi'Hk2] 
+ Fi'\k,)Fi'Hk2)Fi'Hk,)+Fi'Hki)Fi'Hh)Fi'Hk2] 



(Al) 



where the 6 permutations guarantee the symmetry of B^ in its arguments. This leads to Eq. (19). We now expand 
the fci integrand in powers of k/ki up to second order. For notational simplicity, we only consider the contribution 
from a single term a and two permutations (2) -H- (3), and drop the subscript a for the moment. This yields 



d^ki 

(27r)3 
d^ki 
(27r)3 

d^ki 



i7'(2)(k^)^(3)(_k^ - k) + (2) ^ (3) 

F(2)(ki) fi - ¥d, + h^yd.d^ F(^\-\i^) + (2) ^ (3) 
/'(2)(k^)^(3)(_ki) - F(3)(-ki)Fa,F(2)(j^^) ^ ii?'(3)(-ki)fc^Fa,a,i?(2)(]^^) ^ ^2) ^ (3) 



(A2) 



where all derivatives are with respect to ki. Wc now perform an integration by parts for the second and last terms. 
The former term (linear in k) cancels with its permutation (2) -H- (3), yielding 



(27r)3 



2i^(2)(k,)^(3)(_j^^) ^ ii?(3)(_j^^)^»pg^5^.^(2)(i^^) ^ \F^:i'^{)^^)k'¥d,djF^^\-\.{) 



(A3) 



Note that we have not used that F^"^^ (ki) — i^^*-' (fci) so far. We now use this fact however to obtain, defining /i == k-ki. 



k'k^d.djFiki) - T2 [(1 - k^^)kiF'{ki) + M'fc?^"(fci)] , 



(A4) 



where we have denoted derivatives with respect to fci with primes. We thus obtain, up to terms of order k"^ /kf (cubic 
terms drop out just like the linear terms did) 



(*MiMi(2)>=i:{/i|i^<«w.-/^ 



d^ki_ 



2Fi'Hk,)Fi'\k,) 



IP ^r 



+ -^^Pi'Hki)[{i-f^')k^F^^'Hki) + f^'k!F^^^Hk,) 



ip ^, 



+ 2¥/^'^^'''^ [i^-f^')kiK^'Hki) + fi'k!F^^'Hk,) 



+ {(123) ^ (231)} + {(123) ^ (312)} 



(A5) 
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This expression is now in the desired separable form, and the fi integral becomes trivial. We now introduce some 
notation, 



C^i(0 ^ 



4i) 



^V^aLv) 



(27r)3 

2fe) _ 1 



ML{k)F^\k)e'''-'' 



'O.L =^J^ dk,klMl{k^)F^\k^)Fi=\k^) 



27r2 



+ (0 ^ (.?) 



dki 



{Pi^Hk,) [2fci^(^)(A:i) + fc?F^(-''Hfci)] + (*) ^ (j)} 



Note that for all spectral moments, a^^"^^^ = a^^^^\ This allows us to write Eq. (24) in the compact form 

{6l{1)SI{2)) = J2 {2eiM-r^ +41.M4^S +2perm.| , 

where the two permutations stand for the cyclic permutations 123 — > 231, 312. 
We now turn to Eq. (37). This correlator is given by 



(A6) 



(A7) 



{5Ul)y.m = ^{SL{l)SU^)) 



2^2 



d^k 
2alJ JM" 



ML{k)e 



«kr 



d^ki 

(2^ 



X,(|ki|)M,(|ki + k|)B^(|k|, |ki|, |ki + k|) , (A8) 



where Ais{k) = Ws{k)Ai{k). This expression is very similar to Eq. (19), the only difference being that the filter 
functions under the ki integral involve Wg rather than W^. Eq. (A7) then straightforwardly translates to 



r 2(23) , 2(23) N 

(^^(l)y.(2)) = J2M eiir)^ + k^l^(r)^ + 2 perm. , 



(A9) 



where 



T2(ii) = 



27r2 



dk,klMl{k,)F^\k,)Fi^Hk,) 



2(ii) 
'Xas 



127r2 



dki \Ms{ki)Fi'Hki) 



2ki(^MsFi^'>)^ +kl(^MsFi^^^^ 



+ (*) ^ U) 



(AlO) 
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